Topologically protected midgap states in complex photonic lattices 
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One of the principal goals in the design of photonic crystals is the engineering of band gaps and 
defect states. Drawing on the concepts of band-structure topology, I here describe the formation 
of exponentially localized, topologically protected midgap states in photonic systems with spatially 
distributed gain and loss. When gain and loss are suitably arranged these states maintain their 
topological protection and then acquire a selectively tunable amplification rate. This finds appli- 
cations in the beam dynamics along a photonic lattice and in the lasing of quasi-one-dimensional 
photonic crystals. 



anism to induce the midgap state in the beam propaga- 
tion through a photonic lattice; the beam can then be 
manipulated via adiabatic pumping of light. In an alter- 
native realization that includes active components, the 
midgap state constitutes a selectively amplified mode in 
a quasi-one-dimensional photonic crystal laser. 

Complex Su-Schriejfer-Heeger model. — The SSH 
model was originally introduced to describe fractional- 
ized charges in polyacetylene; an exponentially localized 
midgap state then forms at a defect in the dimerization 
pattern [iJl- I consider a version (the cSSH model, 
shown in Fig. [T]) which applies to photonic lattices and 
crystals and incorporates distributed loss and gain [isl - 
17[. The original SSH model consists of a tight-binding 
chain with alternating coupling constants ta and th 
(for being specific let us assume ta > th > 0), and a 
defect in this sequence which supports the topologically 
protected midgap state (see Fig. [1]). The fundamental 
unit cell is composed of two sites (labeled A and B) 
with amplitudes ipr^^ and ipr^\ where the integer n 
enumerates the unit cells. To explore the effects of loss 
and gain in photonic realizations I consider a staggered 
complex onsitc potential ijA = ij + on the A sites 
and i^B ~ ~ ij on the B sites. This modification 
defines the cSSH model. The tight-binding equations 
read 



Since the inception of the field |2| , the design of pho- 
tonic crystals with band gaps and defect states has been 
facilitated by drawing analogies to condensed matter sys- 
tems. A novel impetus for such endeavors is provided 
by the discovery of topological insulators and supercon- 
ductors, systems which occur in distinct configurations 
that cannot be connected without closing a gap in the 
band structure and consequently display robust surface 
and interface states [1, 01 • Recent works have started to 
transfer concepts of band-structure topology to the pho- 
tonic setting. Thus far, this has opened up avenues for 
unidirectional transport 0, Q , adiabatic purnpinff of light 
, and creation of photonic Landau levels d, U as well 
as the creation of bound and edge states via dynamic 
modulation in the time domain IQ, 11 1 . 

The practical utility of topological concepts in photon- 
ics will depend much on the robustness versus absorption 
and amplification. These processes do not have an elec- 
tronic counterpart; they render the effective Hamiltonian 
non-hermitian, and break time-reversal symmetry — but 
in a different way than a magnetic field, whose presence 
of absence enters the topological characterization of elec- 



tronic band structures jl2l Il3j . One may therefore won- 



der whether topological protection can survive the pres- 
ence of gain and loss. 

Remarkably, as shown here for a c omp lex version of 
the Su-Schrieffer-Heeger (SSH) model[l4|, such robust- 
ness can be demonstrated for a photonic realization of 
topologically protected midgap states, localized at an in- 
terface in the interior of the system. Under the infiu- 
ence of spatially distributed gain and loss 15-l3|, these 
states not only maintain their topological characteristics 
but also acquire desirable properties that do not have 
an electronic analogue — the midgap states can be selec- 
tively amplified without affecting the extended states in 
the system. 

The selective amplification of the midgap state can be 
utilized in beam manipulation and lasing. For instance, 
the considered model can be realized as the coupled- 
model theory of a photonic lattice with alternating lat- 
tice spacings. In a setup with passive and lossy compo- 
nents, the midgap state can be rendered lossless while all 
other states suffer identical losses. This provides a mech- 
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where i„ is the intradimer coupling and is the intcr- 
dimer coupling. The infinitely periodic system exists in 
two configurations — a configuration a where i„ = ta and 
t'^ = tb, and a configuration /3 where the values are in- 
terchanged such that tn = tb and t'^ = ta- These config- 
urations are associated with Bloch Hamiltonians 
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delivering identical dispersion relations 



e±{k) = 17 ± 
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FIG. 1. (a) Complex Su-SchriefFer-Heeger (cSSH) chain with ahernating couplings ta and ti, as well as alternating imaginary 
onsite potential i'jA = i(7 + 7) and 173 — 1(7 — 7) (describing loss or gain in the photonic applications). For 71 < the system is 
in the a configuration while for 71 > it is in the /3 configuration, (b) Dispersion Ree(fc) of the extended states, for tb = 0.6 fa 
and 7 = O.Sia. These states have Ime(fc) = 7. (c) Dispersion in the complex eigenvalue plane, including the midgap state at 
eo ~ i^A, which forms due to the coupling defect, (d) As in the original SSH model, the midgap state is exponentially localized 
around the interface and is confined to the A sublattice. (e) Topological characterization of the cSSH model. The outer curve 
shows the trace of the pseudospin vector S(fc) in the xy plane, for the extended states in the upper band e+(fc) of an infinitely 
periodic chain in the a or /3 configuration {tt — 0.6 ta and ^/ta ~ 0.3, 0, or —0.3). As in the original SSH model, Sz ~ 0, so 
that the configurations can be characterized in terms of their winding number (0 in a and 1 in /?). The enclosed loops show 
the trace of the function g{k) whose position in the complex plane determines the direction of S(fc) according to Eq. (O. These 
loops encircle the origin in jS but not in a. 



for extended states with dimensionless wavenumber k. 

In the original SSH model with 7 = 7 = this re- 
sults in two bands, symmetrically arranged about e = 
and separated by a gap A = 2{ta — In the cSSH 
model these bands are shifted into the complex plane, 
corresponding to decaying states if Im e < and ampli- 
fied states if Ime > 0. However, this shift is uniform if 
I7I < 7c = A/2, which is imposed henceforward. Under 
this condition, all extended states experience the same 
overall gain (7 > 0) or loss (7 < 0). In the particular 
case 7 = of balanced loss and gain, the dispersion re- 
mains real, which can be explained by the VT symmetry 
ax[l-L{k)]*a^ = H(fc) with Pauli matrix cr^ 

The midgap state appears when the two configurations 
are coupled together. In Fig. [Ija), the system is in the 
a configuration for n < and in the (5 configuration for 
n > 0, which results in a coupling defect in the middle 
of the sample. The spectrum consists of extended states 
within the two bands, plus an additional state at = 
i'^A, as shown in Fig. [Ub,c). Going back to Eqs. ([T|), 
this value admits an exponentially localized solution with 
= (_i^/i^)-!n| andV-f = [Fig.md)]. In the original 
SSH model the midgap state sits at eo = and preserves 
the symmetry of the spectrum. In the cSSH model the 
midgap state breaks this symmetry in a way that directly 
impacts on its amplification or decay rate — the midgap 
state is more stable than the extended states if 7 > 0, and 
less stable if 7 < 0. This has a topological origin, which 
is discussed at the end of this work. Prior to this I discuss 
applications of the selective amplification mechanism for 
the manipulation of beams and lasing. 



Beam dynamics. — Let us first consider the manifesta- 
tion of the midgap state in the beam propagation along 
a photonic lattice, composed of single-mode waveguides 
as shown in Fig. [5] Experimentally, such lattices can be 
realized using optical fibers, quantum wells, or femtosec- 
ond laser-writing techniques, producing in all cases arrays 
of waveguides with a fixed cross-sectional geometry per- 
pendicular to the propagation direction z [iB-ii]. In this 
setting the parameters 7a and 7s describe the intrinsic 
propagation constants of the waveguides, which are lossy 
if 1A,B < and amplifying if 7a, s > 0. The couplings 
take the values ta and tf,, depending on whether the spac- 
ing between the waveguides is a or 6, respectively, and the 
midgap state now arises from a defect in an alternating 
spacing sequence. Modes with Ime > exponentially in- 
crease along the propagation direction z while those with 
Ime < decay. 

I now set 7A = and 73 = — 27 < 0, corresponding a 
setup with passive A sites and lossy B sites. The midgap 
state is then lossless (eo = 0) while the extended states 
decay uniformly according to Ime = 7 = —7 < 0. Figure 

illustrates the beam propagation in a lattice of 101 
fibers and a spacing defect in the center of the system. 
Panel (a) depicts the arrangement of the fibers close to 
the center of the sample. In panel (b), a broad wave 
packet is fed into the lattice with t?, = 0.2 ia, 7 = 0.05 ta. 
After a short transient the midgap state is populated 
and propagates without attenuation. In panel (c), the 
light is fed into a single A fiber close to the center of the 
sample. Again, the midgap state is populated; it is now 
less localized because here tf, = 0.6ta. In panel (d), the 
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FIG. 2. (a) Realization of the cSSH model in a photonic 
lattice of single-mode waveguides with intrinsic propagation 
constants 7a and 7s as well as alternating spacings a and 
6, and a defect in that spacing sequence (around a; = 0). 
(b) Beam propagation of an initially broad wave packet in 
a lattice of 101 waveguides with tb = 0.2 fa, ')a = 0, 7b = 
—0.1 fa. (c,d) Beam propagation with light fed into an A or B 
fiber close to x = 0, for a lattice with — 0.6 fa, 7a ~ 0, 7s — 
—0.1 fa. (e) Adiabatic pumping of light: Waveguide geometry 
close to the center of the system, (f) Beam propagation in a 
lattice of 101 waveguides, with tt — 0.2 fa, 7a = 0, 73 = 
-0.2 fa. M 



light is fed into a neighboring B fiber of the same lattice. 
The beam quickly subsides as the midgap state is not 
populated. 

Figure [Ue,f ) demonstrates the feasibility of adiabatic 
light pumping 0, I^Si Hi in a lattice where the interface 
gradually shifts by 5 unit cells to the right. In the tran- 
sient region the couplings f„ and fj^ interpolate linearly 
between ta and ff,, with = 0.2 fa, 7 = 0.1 fa- Note 
that the shift of the beam is opposite to the shift of the 
interface. 

These results generalize to systems with 7^ 7^ 0. At 
fixed 7, this situation differs from the passive realization 
by a z-dcpcndent intensity scaling exjp{2"fAz). In active 
realizations with 7^1 = 7 + 7 > > 75 = 7 — 7, I7I < 7c, 
the midgap state is the only amplified state while the 
extended states all decay. 

Laser applications. — When a system with active com- 
ponents is confined in the z direction, the midgap state 
serves as a selectively amplified lasing mode. Figure [3] 
(a) illustrates how such a system could be realized using 
an arrangement of amplifying and absorbing (or passive) 
regions separated by gaps of alternating length. This 
provides a topological realization of microlasing with dis- 
tributed gain and loss [H-iii. While the underlying wave 
equation is second order in time or frequency, in standard 
slowly- varying envelope approximation the eigenvalues of 
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FIG. 3. (a) Realization of the cSSH model in a quasi-one- 
dimensional photonic laser with a staggered arrangement of 
active (A) and lossy (or passive) components (B) in a unit cell 
of size ao. (b) Spectrum of a finite system with 101 modes 
{tt — 0.6 fa, 7 = 0.1 fa), in increasing order of Ree; (i = 
0,1,2, ...,100). (c) Same for complex 7 = (0.1 -f 0.2i) fa 
(cSSH' model). The midgap state remains the most amplified 
state, (d) Implementation in a dielectric medium with regions 
of refractive index n = 1 (passive), ha = 2 — 0.01 i (gain) and 
7iB ~ 2 -f 0.01 i (loss). Selectively amplified midgap states are 
predominantly localized in the gain medium, (e) The depicted 
midgap state (index i = 308) is situated between bands 8 and 
9 of a system of length 40 ao • [sj] 



the cSSH model can now be interpreted as the mode fre- 
quencies LJi = Ei + Q around a large central frequency 

n. 

Panel (b) shows the spectrum of a finite cSSH chain 
with 101 lattice points. For 7 > the midgap state is 
selectively amplified in the time domain and thus will 
win the mode competition. With increasing 7 the las- 
ing threshold then occurs at 7 = —7, at which point 
7^ = 7 + 7 = so that Wo crosses into the upper half 
of the complex plane. For 7 < 0, on the other hand, the 
extended states win the mode competition and become 
lasing at 7 = 0. 

Amplifying and absorbing regions with matching char- 
acteristics pose an experimental challenge. In panel (c) 
this is taken into account via an additional alternating 
real part of the onsitc potential, which is equivalent to 
setting 7 to a complex value (cSSH' model). This shifts 
the real part of the midgap state's frequency but docs 
not affect the imaginary part, which still exceeds that 
of the extended states (the latter now acquire a mode 
dependence). 

Panels (d,e) test the applicability of these predictions 
for an implementation of the laser in a dielectric medium 
with refractive index ua = 2 — O.Olz in the amplifying 
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parts and ns = 2+0. Oli in the absorbing parts of the sys- 
tem. These regions have lengths ao/3 and are separated 
by gaps (refractive index n = 1) of ahernating size ao/12 
and ao/4, where Og is the length of the unit cell. The 
results apply to a system of length 40 oq. Midgap states 
form between the lowest-lying bands that approximate 
the continuum limit. In panels (d,e) this is illustrated 
for the example of bands 8 and 9; the midgap state is 
localized in the amplifying regions and its frequency lies 
much higher up in the complex plane than those of the 
extended states. The results correspond well to the pre- 
dictions of the cSSH model, with minor deviations mostly 
in line with the cSSH' model. 

Topological characterization. — Finally let us discuss 
how the particular features of the midgap state relate 
to the topological properties of the cSSH model. As in 
the SSH model, the difference between the a and (3 con- 
figuration is captured via a topological phase associated 
with the Bloch functions 0, |j, [2^. To formulate this 
characterization it is convenient to write the eigenvectors 
of Hamiltonian ([2]) as 



ip{k) = N 



fi-k) 
e{k) - i^A 



(4) 



where N is the normalization constant. Each extended 
state can then be associated with a pseudospin vector 



2)) — {Sx, Sy,Sz)- 



(5) 



^ , thus 5^ = so that S is 



In the SSH model, \ipA\ = \'Pb 
confined to the a;?/- plane. Remarkably, as long as I7I < 7c 
this property remains preserved in the cSSH model, for 
any value of 7. From the expression given above, the 
direction of the pseudospin can be read off from 



Sx + iSy oc (£(fc) - i^A)f{k) = g{k). 



(6) 



Figure [TJe) uses this relation to trace out the pseu- 
dospin as k passes through the Brillouin zone (— tt < fc < 
TT, shown for the upper band e+(fc)). In the a configura- 
tion, the function g{k) (enclosed loop) does not encircle 
the origin of the complex plane; the pseudospin therefore 
librates and traces out an arc (winding number 0, topo- 
logical phase 0). In the (3 configuration g{k) encircles the 
origin; the pseudospin therefore rotates and traces out a 
circle (winding number 1, topological phase tt). 

This topological robustness is further illuminated by 
considering the chiral symmetry azH*az = —H, where 
az acts in the AB subspace. One can verify directly from 
Eq. ([T]) that this symmetry persists in the cSSH model. 
For the pure configurations, (Tz['H{—k)]*az = — 'H(fc), 
which implies particle-hole symmetry of the dispersion, 
while in the presence of the defect the chiral symmetry 
guarantees the existence of a midgap state [l^l ■ Due to 
the localization on the A sublattice, this state possesses 
a fully polarized pseudospin S = (0,0,1) and inherits 



the complex potential on this sublattice, which thus de- 
termines its eigenvalue eg = i"fA = i"f + i"f. Notably, 
this is still consistent with the constraint Eq = —Eq im- 
plied by chirality; the specific value Eq = only follows 
in the hermitian limit of the SSH model. The extended 
states populate both sublattices equally, which results in 
Ime(fc) = 17. Therefore, the midgap state is more stable 
than the extended states if 7 > 0, and less stable if 7 < 0. 

Conclusions. — In conclusion, photonic systems can 
exhibit exponentially localized, topologically protected 
midgap states whose stability may be controlled via dis- 
tributed loss and gain. Such states can be induced in 
beam propagation through photonic lattices, where they 
provide a platform for adiabatic pumping of light, and in 
photonic crystal lasers with inhomogeneous gain, where 
they exhibit selective level amplification. Remarkably, 
the midgap states maintain their topological protection 
even though the loss and gain renders the underlying 
Hamiltonian nonhcrmitian and breaks the time reversal 
symmetry of the system. This demonstrates the utility 
of topological concepts in genuinely photonic settings. 
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midgap state now arises from a defect in an alternating 
spacing sequence. In this setting the parameters 7a and 
7b describe the intrinsic propagation constants of the 
waveguides, which are lossy if 7a,s < and amplifying 
if 7a, s > 0. The propagation of monochromatic waves 
along z is then described by the coupled-mode equations 
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For systems that are homogeneous in z, the eigen- 
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the tight-binding model ([T]). The results in Fig. [2] fol- 
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Nicholson method. 
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